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Abstract 

Solvent-mediated solid-solid phase transformations often result in the 
formation of a porous medium, which may be stable on long time scales 
or undergo ripening and consolidation. We have studied replacement pro- 
cesses in the KBr-KCl-EbO system using both in situ and ex situ exper- 
iments. The replacement of a KBr crystal by a K(Br,Cl) solid solution 
in the presence of an aqueous solution is facilitated by the generation 
of a surprisingly stable, highly anisotropic and connected pore structure 
that pervades the product phase. This pore structure ensures efficient 
solute transport from the bulk solution to the reacting KBr and K(Br,Cl) 
surfaces. The compositional profile of the K(Br,Cl) solid solution ex- 
hibits striking discontinuities across disc-like cavities in the product phase. 
Similar transformation mechanisms are probably important in controlling 
phase transformation processes and rates in a variety of natural and man- 
made systems. 



1 Introduction 

The re-equilibration of solids in the presence of fluids is of great importance to 
the evolution of both natural and synthetic materials. In most cases, the fluid 
is a liquid solvent, but vapor phase transport in air or other gases may also 
mediate solid-solid phase transformation (Jiang et at 2010). Since Cardew & 
Davey (1985) presented the first analysis of the kinetics of "solvent-mediated 
phase transformations" , it has become clear that such transformations play a 
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key role in the mineralogical evolution of the earth's crust, as well as in a wide 
range of industrial processes including control of cement hydration (Tzschich- 
holz et al. 1996) and production of nanostructured devices (Sun & Xia 2004; 
Xia et al. 2009). Pharmaceuticals and explosives often exhibit polymorphism, 
and they may transform from one crystal structure to another during storage 
and/or processing (Rodrigucz-Hornedo & Murphy 1999, Dharmayat et al. 2008, 
Heinz et al. 2009, Achuthan & Jose 1990, Halcblian & McCronc 1969). This 
is important because the solubility, rate of dissolution, bioavailability, melting 
point and chemical stability of pharmaceuticals depend on the physical form 
and crystal structure of the pharmaceutical, and because the processing be- 
havior, including compressibility and granular flow, varies from polymorph to 
polymorph. Similarly, the performance and stability of explosives varies from 
one polymorph to another. In most cases, crystallization from a melt or solution 
leads to growth of the least stable polymorph (Ostwald 1897, Threlfall 2003), 
which may transform into another polymorph in the presence of a solvent. 

Polymorphism in calcium carbonate is well known because of its impor- 
tance in geological systems and biogeochemistry, and it provides a good example 
of polymorphism and solvent mediated polymorphic transitions in an aqueous 
mineral system. The thermodynamics, kinetics and chemical mechanisms of 
solvent mediated transformations between amorphous CaCC>3, ikaite, vaterite, 
aragonite and calcite have been extensively investigated (Bischoff 1968, Kiyoshi 
Sawada 1997, Mitsutaka Kitamura 2002, Jun Kawano et al. 2009, for example). 
The morphologies of the new mineral phases formed via nucleation and growth 
have been investigated, but this work has been focused on the early stages of 
growth (small scale morphology) and the effects of organic molecules, in order 
to better understand the behavior of CaCC>3 in biological systems. 

Through a variety of recent studies, the mechanism of "coupled dissolution- 
precipitation" (Wang et al. 1995; Nahon & Merino 1997, Putnis 2002) or 
"microphase-assisted autocatalysis" (Anderson et al. 1998) has been verified 
as a key process during mineralogical transformations under a wide range of 
geological conditions (Putnis 2002). Precipitation of a solid product, which is 
more stable than the initial phase, removes dissolved reaction products, and 
this facilitates the transformation by maintaining the solution in an undersat- 
urated state with respect to the initial solid phase. Although most studies 
have focused on the micrometer to nanometer scale characteristics of the in- 
dividual interfaces at which phase transformations take place, there is recent 
evidence that the porosity generation associated with such processes may also 
be linked to fluid migration and advective solute transport. Hence, the mech- 
anisms of solvent-mediated phase transformations may have implications for 
geological alteration and metamorphic processes on very large (up to several 
kilometers) scales (Plumper & Putnis 2009). Regional-scale alteration of rocks 
is often associated with significant mobilization, transport and concentration of 
elements that otherwise only occur at low abundances, thus solvent mediated 
phase transformation is ultimately linked with the formation of economic ore 
deposits (Oliver et al. 1994). 

Very often, the new solid phase grows epitaxially on the old one (Cesare & 



2 



Grobety 1995, Putnis & Putnis 2007), and the replacement is pseudomorphic 
(the replacement product crystal has the habit of the original crystal instead 
of its own usual habit). Although dissolution at one face of the original crystal 
and nucleation on another face would be a plausible replacement mechanism, 
pseudomorphism is more naturally explained by the propagation of a reaction 
front through the original crystal, with dissolution and precipitation occurring 
simultaneously at the reaction front. Porosity is essential for the advancement 
of a mineral replacement reaction front. Without porosity, the reacting phase 
would quickly become "protected" from the fluid phase by a layer of imperme- 
able solid reaction product (s), and the dissolution-reprecipitation process would 
cease. The protection of the original crystal after a thin protective layer has 
formed is similar to passivation during corrosion. In fact, in some corrosion 
processes, the corrosion film is porous, and this facilitates rapid corrosion and 
formation of a thick corrosion product film (Cox & Yamaguchi 1994). In an- 
odization, a thick, highly structured porous film, which is formed on a metallic 
substrate via anodic oxidation (Crouse et al. 2000, Gong et al. 2001), grows 
via dissolution and oxidation at the base of the porous film (Wu et al. 2007). 
Under these circumstances, corrosion/anodization becomes a solid-solid replace- 
ment process with the formation of a porous reaction product. Access of flu- 
ids to reactive surfaces area may be maintained by reaction-driven fracturing, 
and/or pore-formation dominated by growth and dissolution processes. Frac- 
turing depends on volume changing reactions generating sufficient stresses to 
cause fracturing, and it may occur due to expansion or shrinkage (Jamtveit et 
al. 2000; Malthe-S0renssen et al. 2006; R0yne et al. 2008). Pore formation 
has been documented during replacement processes in many systems, including 
alkali feldspars (Worden et al. 1990; David et al. 1995), KBr-KCl-H 2 (Put- 
nis & Mezger 2004; Putnis et al. 2005) and fluorapatite (Harlov et al. 2005). 
Porosity generation by replacement reactions in solution has also been used to 
produce nanostructured metals for a variety of applications, including catalysis, 
optics, electronics and biological sensing (e.g. Sun & Xia 2004). In many cases, 
the porous replacement product is a transient morphology. The large interfacial 
energy associated with the porous structure may drive coarsening and porosity 
reduction. As a result, most evidence of porosity may be erased. However, some 
residual porosity such as fluid inclusions usually remains. 

Despite the geological, biological and industrial importance of porosity gen- 
eration by solvent-mediated phase transformation, the mechanisms generating 
the often highly complex porous structures are still poorly understood. Here we 
present an experimental study of KBr replacement by KC1 in the presence of an 
aqueous solution. Our approach allows direct imaging of the reacting interface 
during porosity formation. We show that the interface develops a highly orga- 
nized structure due to strong coupling between dissolution, precipitation and 
local transport in the fluid phase. The porosity network that forms during the 
replacement process is highly anisotropic and it enables transport between the 
bulk solution and the evolving reactive surface. 
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2 Experimental 



We conducted an experimental investigation of the replacement of KBr crystals 
by K(Br,Cl) solid solutions in the KB1-KCI-H2O system. Initially the solution 
was saturated with respect to KC1 only. The composition of the solid solution 
is described by the formula K[Brx ffBr Cli_jf A , Br ], where Xkbt is the molar 
proportion of KBr in the solid. The kinetics and thermodynamic properties 
of this system are well known (Glynn et at 1990; Pollok 2004). Fig. Q] shows 
the equilibrium relationship between the composition of the solid solution and 
the composition of the aqueous solution (modified after Pollok (2004)). The 
final state of the solid depends on the initial volumes of the solution and of the 
crystal. Replacement of KBr by K(Br,Cl) leads to a reduction of the total solid 
volume because 1) at equilibrium, the solubility of K(Br,Cl) is always greater 
than the initial solution solubility and 2) the molar volume of a K[BrxCli_x] 
crystal is smaller than that of a KBr crystal (Dejewska 1999) with the same 
amount of potassium for all X (0 < X < 1). 

The experimental setup allows direct imaging of the interface and its evo- 
lution during the reaction. One side of the crystals was glued to an overlying 
glass plate with a UV curable adhesive, which prevented reaction from that side 
and provided a transparent path to the replacement fronts/interfaces through 
the glass plate and the KBr crystal (see sketch of the in situ optical views of the 
reaction fronts in Figs. WPD). Because the original KBr crystals were optically 
transparent and isotropic, only the interfaces created optical signals, due to a 
change in refractive index. The advancing reaction front was continuously im- 
aged in situ. The evolution of some crystals was stopped in mid-reaction, and 
these crystals were cut and polished, and images were obtained using ex situ op- 
tical microscopy, scanning electron microscopy (SEM), back-scattered electron 
(BSE) microscopy and quantitative electron microprobe analysis. 

2.1 In situ experiments 

KBr crystals (typical dimensions 2x2x2 mm 3 ) were freshly cleaved from initially 
transparent crystals and immersed in a saturated solution of K + C1~ with a 
volume of V=0.54 mL leading to reaction. The saturated KC1 solution was 
obtained by dissolving KC1 powder in deionized water. The reaction front was 
continuously imaged using Olympus BX-60 microscope with 4X-60X objectives 
and a Cascade CCD camera (spatial resolution 653x492 pix 2 , intensity resolution 
14 bits). White light was used for most imaging and green light (\ g = 547 nm) 
was used for interference pattern quantification. The cavity thickness difference 
between two fringes is A 9 /(4 • n), where n = 1.36 is the refractive index of the 
solution. Different magnifications were used, ranging from 0.34 to 5.1 yum/pixel. 

2.2 Ex situ measurements 

The crystals were removed from the solution, ground roughly under an ethanol 
solution (a few 100 microns of material was removed to access the inner struc- 
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Figure 1: Solubility phase diagram of the KBr-KCl-H^O system (modified after 
Pollok (2004)). The solidus (s eq vs Xkbt) and solutus (s eq vs XBr.aq) curves 
are shown, where s eq is the solubility of the solution and Xs r ,aq is the molar 
proportion of K + Br~ in the solution. The compositions of the different layers 
(XffBrJ X^ Br ) are indicated on the solidus curve. Assuming equilibrium 
conditions, the local solution compositions, shown in the figure (on the solutus 
curve), are deduced. Insert: Expansion of the solubility diagram in the region 
most relevant to the experiments. 

ture) and immersed in acetone in order to remove all traces of water (these 
steps were performed within 2 min). Ethanol and acetone do not dissolve the 
crystals and they dilute the brine solution sufficiently to ensure that no further 
precipitation or dissolution can occur. After 12 hours, the crystals were removed 
from the acetone bath, polished and prepared for analysis. The interface was 
observed from the side (view 2/xz-plane in Fig. [2]) by optical and BSE imaging 
and quantitative electron microprobe analysis. 
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3 Results 



The evolution of the replacement front may be divided into an initial transient 
stage, which lasts approximately 30 minutes, and a "steady state" in which 
the shapes and the compositions of all structures at the replacement interface 
remain essentially invariant while the KBr crystal is examined. These recrystal- 
lization structures have not been observed before because surface energy driven 
recrystallization (without chemical change) erases the complex small scale pore 
structure and leaves the pore geometry observed ex situ by Putnis & Mezger 
(2004) when the crystal is removed from the solution bath. Our approach cir- 
cumvents this difficulty by in situ observations or by freezing the system before 
further ex situ analysis. We first describe the striking "steady state" geometry 
and then the initial transient state that casts light on how the recrystallization 
interface selects the observed steady state. 

3.1 Qualitative optical observations at steady state 

Figs. [2] and [3] show schematic representations of the experimental system and 
several images that illustrate its qualitative behavior. 

The pictures show that the interface between the unaltered KBr crystal and 
the porous K(Br,Cl) saturated with KBr/KCl/H20 solution has a complex ge- 
ometry This interface is bounded by two planar replacement fronts, Ii and 
I<i. I\ separates the region that contains only KBr from the region that also 
contains K(Br,Cl) and solution, and I2 separates the region that contains only 
K(Br,Cl) and solution from the region that also contains some unaltered KBr 
(Fig. [2]ACD). The region between planes 1\ and I2 is the "leading replacement 
zone" in which KBr is replaced with K(Br,Cl). In this zone, replacement occurs 
inside distinctive, generally unconnected, more-or-less cylindrical structures that 
extend from front I\ to front I2, and often some distance beyond I2 (Fig. [2J4) . 
Within the leading replacement zone, the material between the cylinders is 
unaltered KBr, which has never been in contact with the solution (as Fig. Q] 
indicates, once KBr has dissolved in KBr/KCl/H20, it cannot be reprecipitated 
as pure KBr). Observed through the unaltered KBr, from the direction perpen- 
dicular to the interfaces (called view 1), the cylinders appear brighter than their 
surroundings, have cross-sections of varying size (10 2 - 10 4 /im 2 ) and occupy ap- 
proximately 60% of the plane in which they lie (Fig. [2]S). All the cylinder tops 
are almost coplanar and correspond to front Ii (the most advanced boundary 
of the replacement zone) . The observation of interference fringes indicates that 
a thin fluid filled cavity lies at the KBr end of every cylinder (Figs. [2]B, |4]4). 
These fringes (colored in Fig.[2]B, black and white in Fig.HH), visible in reflected 
light due to reflections on both surfaces of the cavities, indicate that the cavities 
become thinner towards their edges (as Fig. [2p also clearly shows). 

Observed in view 1 the disc-like cavities move towards the imaging objective 
(into the KBr crystal) while retaining their gross shapes. This translation of the 
cavities can only occur by dissolution of KBr crystal at their advancing bound- 
aries and precipitation of approximately the same volume of K(Br,Cl) at the 
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opposing interface, which grows the replacement crystal into the advancing disc- 
like cavity. In the center of each disc there is a dark region (Figs. [2]S, |4]A, [5]A) 
with edges that arc aligned with the low index (cubic) crystal planes. The cylin- 
drical structures contain newly precipitated crystals organized into several layers 
separated by brine, which can be seen in views perpendicular to the replace- 
ment interface (view 2, Figs. In transmitted light (view 1, Fig. [5]B) 
the central regions of the cylinders appear brightest, indicating that a fairly 
straight channel or pore (C) connects the disc-like cavities between the newly 
precipitated K(Br,Cl) and extends towards the bulk solution. 

Up to 6 layers of K(Br,Cl) per cylinder structure, which extend above and 
below front I2, were observed (3-4 are visible in the optical images, 6 in the 
BSE image), and they are separated by disk shaped fluid filled cavities. The 
first cavity (rightmost vertical dashed lines in Figs. ^OD,. [6]A) separates the 
original crystal from the first layer of the newly formed crystal, and it is the 
one that is observed in view 1 (Fig. [5]B). All the cavities inside one cylinder are 
similar in lateral and vertical extent to the first one, and they are organized 
around the same channel. 

The channels (C) can be observed in the ex situ image (view 2) when they 
are cut by the slice plane (Fig.[5]C). They are indicated in the sketchs Figs. 04, [31 

The replacement cylinder is always oriented parallel to one of the high index 
crystal planes of the original crystal. Thus, when the original KBr crystal is 
miscut with respect to the high index plane with an angle 7 the angle between 
the plane of the individual cavities (or the K(Br,CL) layers) and the front I2 
(or the original surface) is also equal to 7, and the angle between the axis of the 
central pore and I2 is tt/2 — 7. 

The pore volume inside each cylinder, created by the solid volume reduction 
associated with the replacement process, has a strongly anisotropic structure 
consisting of a central, essentially linear, pore along the axis of each cylinder, 
which connects an array of disk-like cavities, which are oriented perpendicular 
to the central pore. The characteristic scales of this surprisingly ordered pore 
volume are on the order of 10 microns or larger, and the pore can be easily 
resolved by an optical microscope. The scale of the porosity in the rest of the 
K(Br,Cl) (outside cylinders) formed by the replacement reaction is typically 
smaller than 1 (im, which is too small to be resolved with the present setup. 
This fine porous material appears in black in Figs. [2]B CD, (5]4 B , mainly because 
it strongly scatters light. The replacement material in the "trailing replacement 
zone" behind advancing interface I2 (region R in Fig. [2]A) can only be observed 
in high resolution ex situ images (Fig. \5jp) . Dark patches are former fluid filled 
pores (the fluid is evacuated during processing) that have a much more random 
morphology than the channels (C) , and it can be seen that traces of the succes- 
sion of replacement layers remain some distance behind interface I2 in region 
R. 

The distance between the two replacement interfaces I\ and I2 increases 
slowly as they advance through the crystal, and the interface seems to translate 
as a whole in time lapse videos (Video SI). 
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3.2 Quantitative measurements at steady state 

In situ time lapse images and ex situ optical images were used to determine a 
number of spatiotemporal relationships (see Figs. 21 H]). 

The growth of the thickness of the replacement material L(t), the distance 
from the bulk solution to front 1\ (the leading boundary of the replacement zone) 
can be represented reasonably well by L(t) = \[T) ■ t, with D ~ 1.7-10~ 10 m 2 /s. 
This effective diffusion coefficient is an order of magnitude smaller than the dif- 
fusion coefficients of simple univalent electrolytes such as KBr and KC1 in water 
at room temperature, and it is consistent with the idea that the replacement 
of KBr by KC1 is limited by the diffusion of solvatcd K + , Br - and Cl _ in the 
tortuous pore space inside the K(Br,Cl). Although the central pore associated 
with each replacement cylinder is almost straight, the cross-section of the chan- 
nel is much smaller than the area of the KBr crystal exposed to the liquid phase 
along the leading front, Ii. Consequently ion diffusion between the bulk KC1 
solution and the KBr is reduced relative to diffusion in a bulk solution. The 
ion diffusing between the bulk KC1 solution and the KBr are diverted into the 
disk-like replacement cavities and the porous replaced crystal between trailing 
front I2 and the bulk KC1 solution. 

Longitudinal composition profiles of the solid phase from the contact with 
the bulk solution {z = —L(t)) to the front I\ (the leading boundary of the re- 
placement zone, z = 0), obtained from electron microprobe analysis, are shown 
in Fig. [6] Measurements performed at different times can be scaled quite well 
onto a common curve by dividing the positions by L{t). The solid composition 
varies from Xkbt = 1 at the reaction front I\ to Xk Br ~ at the contact with 
the bulk solution. A jump in composition from Xkbt — 1 to 0.72 occurs across 
the first solution filled disc-like replacement cavity. A detailed composition pro- 
file over one replacement cylinder showed that each replacement layer within a 
cylinder has a constant composition. These compositions are denoted by Xj (Br , 
X\ Br , ... (from bulk KBr to bulk solution), and X\ Br = 1. The composition 
jumps from one replacement layer to the next across successive replacement cav- 
ities, and the largest jump is across the cavity adjacent to the original crystal. 
Each replacement layer moves by dissolution on one side and growth on the 
other, but their compositions do not change (Fig. |6]D). This observation is im- 
portant because of the insight that it provides into the replacement mechanism. 
Dissolution and precipitation are coupled by transport of K + , Br~ and Cl~ 
across the thin solution filled replacement cavities, but they do not occur in the 
same place. Due to the fluid composition, dissolution occurs on the rightmost 
part (Figs. [6]AC , [2]CI?) of the replacement cavity, and crystal growth occurs 
on the leftmost boundary. This can be seen on time-lapse images (Video S2): 
Dissolution on one side of the replacement cavity is coupled with precipitation 
on the other side of the cavity via diffusion of ions across the cavity, and the 
coupling of precipitation and dissolution, across the fluid filled cavities, explains 
why the replacement cylinders appear to translate. Each compositional layer 
moves by precipitation on one side and dissolution on the opposite side. To 
maintain a more-or-less constant composition in the solid replacement layers as 



they advance towards the KBr crystal Br~ ions must diffuse through the central 
pores towards the bulk KC1, and chloride ions must diffuse through the central 
pores towards the bulk KBr. 

Seen from view 1, each cavity appears to have a roughly circular shape 
with a radius of r max = 5 — 70/im, which is evaluated from r max = y/A/ir, 
where A is the area of the cavity as seen in view 1. The depths, h(r), of the 
first cavities (those closest to the KBr) can be evaluated using the interference 
fringes observed on Fig. |4]A. The depths are related to the distance r from the 
centers of the cavities by h(r) = a(r max — r), where a = 0.0685 ± 0.0066 is the 
cavity aspect ratio (Fig. HJ. 

By following one cavity for 90 min, we observed that the central channel 
position is essentially fixed with time (Fig. 0J). This illustrates the surprising 
stability of the shapes generated by the coupled growth and dissolution process. 

3.3 Transient stage 

The cylinders, made of solid layers, each with a hole in its center, separated by 
thin fluid filled cavities, are formed during an initial transient stage. The major 
changes (the complex dynamics leading to the formation of the channels and the 
solid layers with different compositions) occur during the first few minutes, and 
this is followed by a slow relaxation to the steady state cylindrical shapes, which 
takes about 30 min. Here we describe the very early part of the transient, which 
is the key to understanding the stable pattern observed in the steady state, based 
on time lapse video images in view 1 with reflected light (Video S3). 

In view 1, light is reflected from both the KBr-brine interface closest to 
the viewer and from the adjacent brinc-K(Br,Cl) interface further away (on the 
opposite side of the fluid filled gap). When these two interfaces are closer than 
5-7 /xm, the two reflections interfere and the interference fringes are used to 
interpret the relative angle and flatness (or waviness) of the two crystal-brine 
interfaces. Dark regions, the absence of reflection, are due to pores or contact 
between the two crystal surfaces on opposite sides of the fluid filled gap. 

Immediately after the brine is added to the KBr crystal small (7-10 /im), 
roughly cubic K(Br,Cl) crystallites grow on its surface. They are all roughly 
aligned and they cover almost half the KBr surface (Fig. W0)- The shape, align- 
ment and lack of reflection from these crystallites suggest that they grow epi- 
taxially on the KBr surface. 

The dissolution front evolves from the space between the crystallites and 
undercuts the contacts between the original KBr crystal and the K(Br.Cl) crys- 
tallites. Meanwhile, the crystallites grow laterally, connect and form larger 
crystals. The precipitation front now becomes a wavy, continuous surface of in- 
terconnected K(Br,Cl) crystals that are detached from the original KBr crystal 
(Fig.ER. 

The dissolution front remains wavy while the original KBr crystal and new 
K(Br,Cl) crystals are out of contact, suggesting that dissolution is transport 
controlled. However, the growing K(Br,Cl) crystals maintain high index surfaces 
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during growth, and the distance between the two crystal surfaces does not exceed 
10 jUm. 

As the pores in the growing K(Br.Cl) crystal become narrower and deeper, 
dissolution and growth slow down and high index surfaces become increasingly 
dominant. 

A few crystals start growing more quickly than the rest, and they keep up 
with the dissolving KBr surface while the other crystals lag behind (Fig. [7]D). 

The more rapidly growing crystals have high index surfaces with misfits of 
up to 3-4 degrees. The distances between the centers of these rapidly growing 
crystals vary from 10 to 100 /jm. The rapidly growing crystals grow both later- 
ally and towards the original, receding, KBr crystal and end up contacting both 
the original KBr crystal and each other (Fig. [7]B) . Open pores are maintained 
at junctions between the fast growing crystals and persistent thin fluid cavities 
develop from these pores (Fig. [7]/ 7 '). 

Approximately 6.5 min after the beginning of the reaction, replacement 
cylinders, each with several layers, are well formed and the layer compositions 
have already reached their asymptotic values (Fig. 

After this stage, there is almost no additional change in the shape of the 
channel, but the shape of the cavity still evolves during the next 10-30 min. 
The discs (immediately behind front Ji) then covers approximately 60% of the 
plane of the original KBr surface (Fig. H]£>). 

4 Discussion 

The results presented above provide important insights into the mechanisms of 
solvent-mediated phase transformation involving solid solutions. The formation 
and maintenance of connected porosity provides transport pathways from the 
bulk solution to the original crystal. This pore-structure prevents drastic re- 
duction of the replacement rate that would otherwise occur due to formation 
of a "protective" layer of product phase, and it maintains a reactive interface 
between the precursor crystal and the solution. The dynamics of the replace- 
ment front quickly evolves towards a quasi-stcady-statc after a transient period, 
which is the key to the entire replacement mechanism and rate. 

During the initial stages of the replacement process, the product phase forms 
7-10 fjLin crystallites by heterogeneous/epitaxial nucleation and growth on the 
reactant (KBr) phase. Although a continuous epitaxial film of the product 
never forms on the reacting surface, this determines the crystal orientation of 
the product phase throughout the entire growth period by subsequent epitaxial 
precipitation on the crystallites, which grow to become the precipitated phase. 
Thus, although most of the newly formed crystal has never been in contact 
with the parent phase (this is true in the transient stage, but also in the steady 
state since fluid- filled cavities lie between parent and new formed crystals), its 
lattice orientation is dictated by that of the parent phase. In this way, the strain 
energy associated with epitaxial growth across an interface of significant lattice 
mismatch remains small. 
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The formation of a porous network, which secures effective transport to 
the reactive surface, is sustained in this system, due to the overall reduction in 
solid volume associated with the replacement reaction (at least ~13% (Putnis & 
Mezger 2004; Pollok 2004)). The early channel formation occurs because some 
crystallites grow faster than others in the brine filled space between the original 
crystal and the replacing crystal (Fig. 0D). This closes the space between re- 
placing crystals except where these crystals merge (Fig. [7]B). This competition 
between the replacing crystals would have to be addressed in detail to under- 
stand the transition from the early stage growth of faceted crystals to the more 
complex growth of porous morphologies in natural systems and elucidate which 
minerals and conditions will result in a replacement cylinder structure similar to 
that displayed by KBr-KCl at room temperature, and which will evolve in other 
ways. After the first disc-like cavity has been formed (Fig. [7]F) "nucleation" 
and development of the second and subsequent cavities could not be directly 
observed. We speculate that the formation of the second and subsequent cav- 
ities is similar to the formation of the first cavity, by nucleation of crystallites 
onto the preceding layer, but confirmation of this hypothesis would require a 
quite different experimental approach. 

In the steady state, there seems to be two main reasons why the pores re- 
main open in roughly the same position: There is a solid volume deficit in the 
replacement of KBr by K(Br,Cl) and the misfit (up to 3-4 degrees) between the 
crystallites that grow to become pore walls hinders complete closing of pores 
and stabilizes their positions since this maintains flat high index surfaces and 
straight pores. Straight pores maximize the rate of solute transport, and the 
rate at which the replacement front advances. The range of the composition 
profiles is expected for a transport-limited reactive system with KBr-rich crys- 
tals in the leading part of the replacement zone and KCl-rich crystals close to 
the bulk solution (and far from the KBr). Indeed, the solution becomes richer 
in Br and poorer in CI when it penetrates the porous material, and assuming 
that the precipitated solid is locally in near equilibrium with the solution, a 
solution richer in Br is expected to precipitate a solid that is richer in Br (from 
the diagram of Fig. [1]). But detailed profiles show very distinct composition lay- 
ers, at least inside the replacement cylinders (the fine-porous material formed 
behind the advancing KBr replacement zone is not resolved with the present 
setup and the variations in composition in the wakes of the cylinders is too 
small to discriminate) , which are due to the discrete nature of the replacement 
that occurs between the two sides of a fluid filled cavity. The composition of 
each layer was shown to be stable over time as the front progressed into the 
crystal, suggesting that epitaxial growth is the main mechanism that occurs on 
the precipitation side of a fluid filled cavity. This suggests that atomic layers 
are deposited onto the replaced crystal with exactly the same composition and 
lattice properties. This prevents a mismatch between the substrate and the de- 
posit that would lead to a strain energy cost. This mechanism is maintained by 
the specific shape of the fluid-filled cavity (with a thickness that is significantly 
smaller than its radius (Fig. [4])), which facilitates the intimate coupling between 
dissolution on one side and precipitation on the other side of the cavity. This 
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mechanism is also very efficient: The composition of the first replacement layer, 
after the first cavity, is KBro.72Clo.28 (28% of the bromine in the KBr has been 
replaced by CI), almost 50% of the Br has been replaced in the fifth replace- 
ment layer, and more than 50% replacement occurs in the sixth replacement 
layer and the porous K(Br,Cl) layer. The different layers form after only a few 
minutes, and assuming that epitaxial growth was efficient at the earliest stages, 
it can be assumed that the layer compositions are determined by the composi- 
tion of the crystallites that were first able to nucleate onto the original crystal 
(or subsequent layers) during the very early stages of the reaction. 

The scaling of the composition-position data onto a common curve (Fig.[S]Z?) 
with the replacement thickness L(t), with L(t) ~ i 1 / 2 , is typical of diffusion 
controlled systems. But once the entire original material has been replaced, the 
replaced material might ripen or equilibrate which would erase or rearrange the 
internal porosity. In most cases, the final product will be a transparent, non 
porous crystal, as observed in the KBr-KCl-H20 system (Putnis et al. 2005), 
or another homogeneous material without a trace of the "transient" porosity. 
Some remnants of the channels might occur as fluid inclusion trails, which are 
often observed in minerals (Crawford & Hollistcr 1986). 

5 Conclusion 

In situ and ex situ studies of replacement processes in the KBr-KCl-H^O sys- 
tem reveal the development of striking pore structures. The fluid film between 
the original and replacement crystals was observed for the first time, and this 
provides the key to understanding the coupled dissolution and reprecipitation 
processes in solvent-mediated phase transformation. Replacement of KBr by a 
K(Br,Cl) solid solution that becomes more Cl-rich towards the bulk solution, 
takes place in an array of discrete more-or-less cylindrical compartments near 
the advancing KBr dissolution front. A complex and highly anisotropic pore- 
structure evolves as the front advances. The rate at which the front advances 
is controlled by diffusion in an anisotropic channelized porosity network, which 
is generated by precipitation of replacement reaction product, with a reduction 
in solid volume. Initially replacement is driven by dissolution and the epitaxial 
growth of replacement product crystals with well defined compositions, influ- 
enced by fluid composition and lattice mismatch effects. At later stages, the 
same mechanisms coupled with solid volume reduction and partial restriction 
of long-range transport allow a surprisingly stable partially ordered pore struc- 
ture to form in a solid matrix with large composition changes that are strongly 
correlated, and mechanistically linked, with the ordered pore structure. The 
replacement preserves the gross morphology of the crystal and the crystalline 
orientation, which are characteristics of such transformations (Putnis & Putnis 
2007), despite the fact that dissolution and reprecipitation occurs with an in- 
termediate liquid film. We have explained how the initial stages of the process 
may imprint the crystal orientation. Our findings have major implications for 
understanding transformation processes in many systems in which reduction in 
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solid volume occurs during epitaxial replacement. Based on the results reported 
here, it is not possible to conclude whether or to what degree the distinctive 
pore structure observed in the KBr-KCl-H^O system is generic. It appears that 
volume reduction and epitaxial growth are necessary for the characteristic be- 
havior observed in the KBr-KCl-H^O. However, a large body of additional work 
will be needed to classify porous solid-solid replacement processes and develop a 
predictive understanding of the evolution of porous morphologies during solvent 
mediated replacement reactions. 
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S: (dark blue) solution 

O: (transparent) original KBr crystal 

R: (light gray) replaced K(Br,Cl) crystal 

Ji: advancing boundary 

of the leading replacement zone 

I 2 : trailing edge 

of the leading replacement zone 

C: pore channel 

D: disk-like cavity 



Figure 2: Optical images of the replacement interface. A) Representation of 
the experimental system with solution (S), replaced crystal (R) and original 
KBr crystal (O). The layered cylinders of organized replacement material ex- 
tend through the leading replacement zone between the front 1\ and the front 
I2 and into the brine saturated porous K(Br,Cl) trailing replacement zone that 
contains no pure KBr (the slices on the edges indicate that the cylinders ex- 
tend into the fully replaced crystal region). If a slice cuts a cylinder along its 
axis a pore channel (C) connecting 1\ to the bulk solution is observed. The 
replacement fronts 1\ and I2 are imaged in two orientations: view 1: image of 
the replacement front planes {I\ and I2), looking towards the x — y-plane, and 
view 2: image perpendicular to the replacement front planes, looking towards 
the x — z-planc (or y — z-planc). B) In situ reflected light microscopy color 
image in view 1 focused on one cylinder top of front I\ . C) Ex situ optical mi- 
croscopy image in view 2 showing the original, transparent crystal to the right, 
the dark, replaced crystal to the left and the two replacement zone boundaries 
I\ (4 cylinders) and li- D) In situ transmitted light microscopy image in view 
2 of a similar region. One replacement^dylinder is observed, with its pore chan- 
nel (C) and at least 3 disk-like cavities (ED). In D) only, the initial crystal was 
miscut leading to a different orientation for I\ and I2 (see text for explanation) . 
Scale bars are 100 /im. 




Figure 3: Unsealed diagram in view 2 of the replacement morphology. Colors 
and symbols refer to the caption of Fig. [2j 
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Figure 4: Characterization of fluid filled cavity morphology. A) Optical image, 
with a green filter, of the first replacement cavity (view 1); the scale bar is 
100 /im. B) Superposition for different times (between t = 30 and 120 min) of 
the edges of the cavity and the central channel (green in A)). C) Reconstruction 
of the thickness profile from optical images over different cavities and times (the 
profile corresponding to picture A) is displayed in red). A linear interpolation 
leads to an average slope value a = 0.0685. This value gives the aspect ratio 
between the thickness and the radius of the cavity. D) statistics of the local 
slopes. The average slope, a, is displayed in dark blue. 
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Figure 5: Pore channel geometry. A) In situ reflected light microscopy image in 
view 1 focused on front I\ . B) In situ reflected and transmitted light microscopy 
image with the same focal plane and region as ^4). Pore channels transmitted 
light directly through the whole replaced crystal. C) Ex situ BSE microscopy 
image in view 2. 4 replacement cylinders followed by 4 pore channels (horizon- 
tally oriented black shapes) can be seen^ Symbols refer to the caption of Fig. [51 
Scale bars are 100 /im. 




Figure 6: Solid solution composition measurements. A) BSE image of a re- 
placement cylinder, scale bar is 100 /im. B) Composition along the axis of the 
cylinder versus distance from first replacement front I\. C) CI concentration 
field from microprobe measurement (intensity proportional to CI concentration). 
The positions of the disk-like cavities separating adjacent layers are indicated 
by the dashed red lines. D) Composition of the different layers inside a replace- 
ment cylinder versus time. E) Composition profile over the whole replacement 
crystal at different times. The whole profile displays composition jumps from 
Xkbt — 0.5 (the composition of the identifiable solid displacement layer near- 
est to the bulk KC1) to Xkbt = 1, which correspond to successive replacement 
layers. I\, z = 0, is chosen as the origin of the z-axis. 
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Figure 7: Transient, early stage, replacement processes. Right: Snapshots of the 
interface (view 1) taken at different times: before immersion (A), 0.2 (B), 5 (C), 
60 (D), 75 (E) and 115 s (F). Scale bar is 100 /im. A close-up of (B) is displayed. 
Left: 2D representations (unsealed) in view 2 of the view 1 pictures. The colors 
of the solution, the replaced and original crystals are identical with those of 
Fig. [2]A. Although the crystals are depicted as separate in the 2D sketch (C) 
and (D) the individual crystallites are connected in the third dimension before 
detaching from the parent crystal. Symbols refer to the caption of Fig. [2] 
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